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Abstract 

We present a theory of electronic properties and the spin blockade phenomena in a gated linear 
triple quantum dot. Quadruple points where four different charge configurations are on resonance, 
particularly involving (1,1,1) configuration, are considered. In the symmetric case, the central dot is 
biased to higher energy and a single electron tunnels through the device when (1,1,1) configuration 
is resonant with (1,0,1), (2, 0,1), (1,0, 2) configurations. The electronic properties of a triple quantum 
dot are described by a Hubbard model containing two orbitals in the two unbiased dots and a 
single orbital in the biased dot. The transport through the triple quantum dot molecule involves 
both singly and doubly occupied configurations and necessitates the description of the (1,1,1) 
configuration beyond the Heisenberg model. Exact eigenstates of the triple quantum dot molecule 
with up to three electrons are used to compute current assuming weak coupling to the leads and 
non-equilibrium occupation of quantum molecule states obtained from the rate equation. The 
intra-molecular relaxation processes due to acoustic phonons and cotunneling with the leads are 
included, and are shown to play a crucial role in the spin blockade effect. We find a quantum 
interference-based spin blockade phenomenon at low source-drain bias and a distinct spin blockade 
due to a trap state at higher bias. We also show that, for an asymmetric quadruple point with 
(0,1,1), (1,1,1,), (0,2,1), (0,1, 2) configurations on resonance, the spin blockade is analogous to the 
spin blockade in a double quantum dot. 
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I. INTRODUCTION 



Gated quantum dots (QDs)^"- with controlled electron numbers are a testbed for probing 
fundamental many-body physics as well as a promising platform for building spintronics 
and quantum information processing (QIP) devices.- Until recently, most experimental and 
theoretical investigations of quantum circuits based on electron spin focused on the single 
and double quantum dot (DQD) devices.-!^. Many essential tasks for operating a qubit have 
been demonstrated in DQDs. For instance, coherent manipulation and readout of one^ and 
two^ spin states have already been experimentally achieved using spin blockade.-'ii'i^ In 
DQDs, spin blockade is used to detect spin using spin-to-charge conversion. For instance, 
the (0, 2) charge configuration cannot be obtained from the (1, 1) configuration if the electron 
spin in the left dot is parallel to the electron spin in the right dot. Detected charge on the 
right dot depends on the relative spin orientations of the two electrons. Thus, spin blockade 
detects spin states (triplet or singlet) of the two electrons in transport spectroscopy or 
charge sensing measurement .-"ii^i^ A physical signature of spin blockade at the triple point, 
(0,1) — 7- (1,1) — (0,2), is the current rectification under different bias directions. In 
positive (forward) bias direction, triplet states will not be populated, and the system does 
not manifest negative differential conductance. In negative (reverse) bias direction, current 
suppression is pronounced once the transitions to the (1, 1) triplet states become accessible 
in the transport window. 

A nontrivial extension of the quantum circuit based on electron spin is the triple quantum 
dot (TQD) with one electron each. This can be appreciated by the comparison of the 
quantum optical properties of a two-level versus three-level systems. Charging and transport 
spectroscopy experiments^ii^"— on the TQDs have already mapped out the stability diagram 
of the devices down to a few electrons. Recent experiments^ii^ have also demonstrated 
coherent manipulations of electron spins in TQDs. The electronic properties of a TQD have 
been investigated theoretically, including topological Hunds rules,— spin-selective Aharonov- 
Bohm oscillations,—'^ the implementation of a coded qubit,—"— voltage-controlled spin 
manipulation,— entangled GHZ state generation,— non-Fermi-liquid behaviour—"— in 
a triangular TQD as well as coherent tunneling adiabatic passage (CTAP) processes for a 
single electron in a linear triple quantum dot (LTQD).— All these theoretical predictions 
as well as quantum information processing in a TQD require an ability to spectroscopically 
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detect spin by, e.g., spin blockade. 

In recent experiments Granger— et. al. and Laird^ et. al. carried out transport 
spectroscopy and charge sensing measurement on a LTQD molecule with one electron in 
each dots. This configuration, denoted by (1,1,1), was tuned to be resonant with the 
two electron configuration (1,0,1). It was assumed that transport proceeded through 
{(2, 0, 1), (1, 1, 1), (1, 0, 2)} resonant configurations, which implied that the central dot was 
biased to higher energy. The presence of doubly occupied dots in the configurations makes 
the Heisenberg model of localized spin configurations inapplicable and a microscopic model 
is required to study the electronic and transport properties of this TQD system. 

Here we extend our earlier theory of a TQD^i^SiSii^ to biased linear molecule at quadru- 
ple points (QPs) and describe spin blockade as a spectroscopic tool allowing the readout 
of electron spin. We analyze the electronic and spin properties of a LTQD as a func- 
tion of energies of each dot within a single-band or multi-band Hubbard model. The 
knowledge of the wave functions of a single-band Hubbard model allows for the qual- 
itative understanding of the low-bias transport through the device, but including more 
than one orbital in the dot will be shown to be crucial for spin blockade. Two differ- 
ent QPs involving the (1, 1, 1) configuration are considered: (a) symmetrical QP (SQP) 
with (1, 0, 1), (2, 0, 1), (1, 1, 1), (1, 0, 2) configurations on resonance, and (b) asymmetrical 
QP (AQP) with (0, 1, 1), (1, 1, 1), (0, 2, 1), (0, 1, 2) configurations on resonance. For SQP, 
the transport goes through (1,0,1) — )■ (2,0,1) — t- (1,1,1) — )• (1,0,2) channels, while 
(0,1,1) (1,1,1) (0,2,1) (0,1,2) is the transport channel for the AQP. Current 
is calculated in sequential tunneling approximation between the TQD and the leads, using 
rate equationsi^i^ to calculate the non-equilibrium steady state occupation of TQD states 
with a source-drain bias. We use Fermi's Golden Rule to calculate the transition rates be- 
tween TQD states by adding or removing an electron due to the coupling between the TQD 
molecule and the leads, and also the transition rate between TQD states with the same 
number of electrons due to the interaction with acoustic phonons.— 

The plan of the paper is as follows. In SecJTl] and Sec JII At we describe the system, the 
Hamiltonian, and the electronic properties of a TQD as a function of detuning A of the 
central dot. In Sec JIIBt our approach to the transport based on the sequential tunneling 
between the leads and the TQD molecule and rate equations are explained in detail. The 
transition rates due to different mechanisms are also discussed. In Sec JIII A] we present 
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results of current calculations for the SQP, and discuss the mechanism of quantum spin 
blockade at low bias. In Sec lIII B| we present results of transport calculations for conven- 
tional spin blockade at the SQP under high source drain bias and at the AQP, and discuss 
how the system at the AQP can behave qualitatively as a double dot around a similar triple 
point with (0, 1), (1, 1), (0, 2) configurations. A brief conclusion is given in Sec JIVI 

II. MODEL 

Figure [1] presents a schematic diagram of a LTQD in contact with the two semi-infinite 
leads and the energy levels of the single QD orbitals. The metallic leads are modelled by one 
dimensional tight binding chains. Each quantum dot, defined by metallic gates on top of 
GaAlAs/GaAs heteroj unction and represented here by a circle contains a controlled number 
of electrons, e.g., one electron each [(1, 1, 1) configuration] in (a) and (1,0, 1) configuration 
in (b). Electrons can tunnel between dots 1 and 2, and between dots 2 and 3, but there is 
no direct tunnel coupling between the two edge dots. Figure [T]^c) shows the single particle 
levels of the individual dots in the LTQD without interdot tunneling. The lowest energy 
bars denote S orbitals (the ground orbitals) in each dot. The energy of the central dot is 
raised by an applied voltage A. This bias can be used, for example, in order to localize 
the two electrons in dots 1 and 3 as shown in Fig. [I](b). In this study. A, comparable 
to Coulomb repulsion U, is used to bring the configurations such as (1,1,1) and (1,0,2) 
on resonance as shown in Fig. [2l We find it is essential to include the excited states, P 
orbitals, in dots 1 and 3 in order to properly account for the transport properties at the 
SQP. The energy separation, A^p, between 5* and P orbitals may also be comparable to A. 
Thus, the electronic properties of a LTQD are described by a multi-band Hubbard model 
with parameters derived from a microscopic Linear Combination of Harmonic Orbitals- 
Configuration Interaction (LCHO-CI) approach for given voltages on the gates. With Cjo- 
(cj^) denoting annihilation (creation) operators for an electron with spin a on orbital i, the 
five-level Hubbard Hamiltonian reads: 

5 5 5 ^5 

Hd= ^ Ei{Vsd)nia + ^ Ujc\^Cj„ + ^ Uihi^fii^ 2 ^ ^ijPiPi^ (1) 

i=l,a i,j = l,cr i=l ij = l 

where EiiVsci) is the sorce-drain bias dependent energy of orbital i, and tij, Ui, and Vij are 
tunnel coupling, on-site and off-site Coulomb repulsion between orbitals i and j respectively, 
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hia = cl^Ci^, and pi = J2a^io-- We assign indices i = 1,2,3 to orbitals of dots 1, 2 and 
3 respectively. The indices i = 4,5 denote excited P orbitals. We will consider only single 
excited orbitals in both dot 1 (i=4) and dot 3 (i=5) for the TQD molecule at SQP. For the 
AQP case with (0, 1, 1) base configuration, the excited orbitals are in dot 2 (i=4) and dot 3 
(i=5). The excited orbital in the biased dot does not play any significant role. 

The TQD device is connected to left and right leads (r = L,R) as shown in Fig. [It^a). 
Electrons in the leads fill up the noninteracting states of semi-infinite tight-binding chains 
with a bulk dispersion relation er{k) = 2tr cos (ka) up to a Fermi level Pl(r), where tr is the 
tunnel coupling between the sites on lead r, a is the distance between sites of tight-binding 
chain, and k denotes the mode of the plane wave for single particle states in the chain. The 
interaction between the leads and the device is modelled as, 

^'■^ = E E (^l(^)4.c... + h.c.) , (2) 

ir,cr k 

A2-Kkamr 



where t\{k) = t[e* '^™'''/v27r is the tunnel coupling between the mode k of the r = L{R) 
lead and orbital v localized in the left dot (r = L) or the right dot (r = R) and in the 
exponent of is 1 for r = L and —1 for r = R. creates an electron with momentum 
k and spin a in the lead r. In this study, tf = for orbitals not in the right edge dot and 
= for orbitals not in the left edge dot. 

Interactions with phonons have already been shown to be important to understand the 
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for in- 



incoherent transport properties of double quantum dots at high bias in Ref. 
stance. We include interaction of electrons in the LTQD with bulk longitudinal acoustic 
(LA) phonons via deformation potential as the mechanism of phonon-induced relaxation at 
low temperature. The electron-phonon interaction Hamiltonian reads, 

5 

He-vh = E E + ^-q) 4c,-, (3) 

where q is the phonon momentum, i and j are TQD orbitals, and ^q(^q) operator annihilates 
(creates) a phonon with momentum q. Mjj(q) = A(g) j '?/'j(r)* exp(— zq ■ r)^/'j(r) is the 
electron-phonon scattering matrix element, ipi{r) is a single particle wave function, and 



^(Q') = Y deformation potential -D, GaAs mass density p, and speed of sound in 

GaAs. The phonon scattering matrix element, Mjj(q), depends on the single particle wave 
function ^i(r) which is obtained from the LCHO^ formalism. 



A. Electronic Properties of a LTQD 



The electronic properties of a triangular TQD molecule with all three dots on resonance 



14. We focus here on the linear 



for = 1 — 6 electrons have been described in detail in Ref . 
molecule where there is no tunneling between the end quantum dots and on the effect of 
detuning A of the energy of the central dot. While in numerical calculations we retain all 
five levels, we retain only the three lowest energy states in this semi-analytical discussion 
of the low energy spectra. For the Hubbard parameters, we set Ui = U, ti2 = ^23 = t, 
ti3 = 0, Vi3 = V, and V12 = V23 = V' . For the on-site energies, we restrict our attention to 
El = E3 = E and E2 = E + A. Since E is just an overall shift in energy, we will simply set 
E = until we explicitly state otherwise. The Hubbard Hamiltonian commutes with total 
S"^ and Sy , so we consider spin-resolved subspaces in the Hilbert space. 

First, we focus on the single particle molecular states of the TQD. We consider the 
Sy = 1/2 subspace and use a localized basis {|1), |2), |3)} , where \i) = 4^|0). In this basis, 
the Hubbard Hamiltonian, Eq.([T]), reads, 

t 

Hie= t A t . (4) 
t 

By inspection we see that a state \D) = (|1) — \3))/\/2, with an energy E^ = 0, is an 
eigenstate. In this state, an electron does not occupy the central dot. This state can block 
the transport in a setting where dots 1 and 3 connected to the source and dot 2 connected 
to the drain, and hence is called a dark state,— in analogy to the coherent population 
trapping in quantum optics. The existence of a dark state can be detected by transport 
spectroscopy^ of an empty dot. As the transport window determined by the applied source- 
drain voltage Vsd is large enough to allow the added electron to enter a dark state, a negative 
differential conductance should be observed in the experiment. Furthermore, Greentree^^*^ 
et. al. proposed to implement CTAP to move an electron from dot one to dot three without 
passing through dot two and for quantum information transfer for a double-dot charge qubit. 

There are two states orthogonal to the dark state \D): the bright state \B) = (|1) + 
\3))/\/2 and the central state \C) = |2) . The 2-by-2 Hamiltonian matrix spanned by 
the bright and central states can be analytically diagonalized, and the two eigenstates are 
expressed as a linear combination of the bright and central state: |Mi) = cos(0)|i?) + 



sin(0)|C) and IM2) = - sm{(f))\B) + cos{(j))\C) , where tan(2(/)) = -V2t/A. We note that 
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Tuning t mostly 



tuning (p allows us to recover Jacobi eigenstates discussed in Ref,. 
controls the amount of mixing between the bright and central state in the two eigenstate |Mi) 
and IM2), whereas tuning A can control the energy spacing between the two eigenstates. 
The energies associated with the three eigenstates \D), |Mi), and IM2), are Ed = 0, Emi = 
(A - At)/2, and Em^ = (A + At)/2 where A^ = VA^TsF. We note that |Mi) is always 
the ground state. 

Next, we address the two-electron case. The simpler case to analyze is the triplet Sy = 1 
subspace, which contains two spins up in the LTQD. There are 3 basis vectors {|Ti), IT2), 
IT3)}, where |Ti) = c^^c|^|0), \T2) 
Hubbard Hamiltonian in this basis reads 



--St 



c|-j.|0), and IT3) = Cg^Cg^ |0), respectively. The 
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A + V t 
t V t 
t A + V' 



(5) 



The triplet Hamiltonian, Eq.(l5]), and the single particle Hamiltonian, Eq.dl]), have the iden- 
tical matrix structure. Therefore, there is a dark triplet eigenstate |Td) = (|Ti) — \T-i))/\/2 
and the bright \Tb) = (|Ti) + \Ts))/\/2 and central state, \Tc) = IT2). Rotating the Hamil- 
tonian, Eq. ([5]) into the basis of bright, central, and dark states, a 2-by-2 Hamiltonian 
matrix coupling the bright and central states is derived. The two eigenstates of the triplet 
subspace are \Mf) = sin(0)|TB) + cos{(j))\Tc) and |Mj) = — cos(0)|Tb) + sin(0)|Tc), where 
tan(20) = \/2t/Ay, and A^ = A + V' — V . The corresponding eigenenergies of the three 
states are 



Exr, 



A + V 



E,,T = A + V'-^(^A, + V(A.)2 + 8t2 
Emj = a + - ^ (a, - V(A.)' + 



(6a) 
(6b) 
(6c) 



The ground state \M'[) is predominantly characterized by \Tc) = \ tits) with spins up in 
dots 1 and 3 because the corresponding coefficient sin(20) ^ 1 — when t/A is small. 
Nevertheless, {M^) still has non-zero presence in both |Ti) = | tit2) and IT3) = | t2t3) 
configurations. In later sections, we will explain how the low bias spin blockade formation 
is related to the small yet finite components of |Ti) and IT3) in \M^) wave function. We 



designate three ground states in each of the spin-resolved triplet subspaces with (5* = 1, Sy = 
1,0, -1) as I T+), I T°), and | T"), respectively. | T+) = | Mf) as was shown above, and 
I T°) is obtained by flipping one spin and performing symmetrization of the wavefunction 
and I T^) is obtained by flipping both spins from | These states will play the major 

roles in the transport through LTQD at the low source-drain bias Furthermore, we find it 
also useful to represent |T+) as: 

\T^) = 1 (\T,) + ^ \T,) + ^ IT,)] , (7) 
V 7 7 / 

where the coefficients |7i(2)| ^ 1- 

Next, we analyze the Sy = singlet state for 2 electrons. We define the follow- 
ing basis 15*2), jS's), 15*4), IS^), \Sq)}. The singly occupied configurations are, \Si) = 

71 (cu4t + ^lAt) |0)> l-^a) = 75 (4^4t + ^lAt) 1*^3) = 75 (clAt + 4;4t) |0)- 

The doubly occupied configurations are, 15*4) = c^cj^ |0), l^s) = 44,4tl'-')' 1*^6) = 
4i4tl0)- '^^^ Hubbard Hamiltonian in this basis reads. 



A + V' 
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2A + U 
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u 



(8) 



The 3-by-3 upper left block, spanned by {15*1), 15*2), I'S's)}, is identical to the triplet Hamilto- 
nian, Eq. ([5]). For |A| small compared to on-site Coulomb repulsion U, the energy spectrum 
of the singlet subspace can be divided into the bands of singly occupied and doubly occu- 
pied configurations, with a gap of the order of U. Under such condition, the energies and 
wavefunctions of the first three lowest singlet states are very similar to those of the triplet 
states, and the mixing between singly and doubly occupied configurations leads to a t — J 
model.— However, if |A| is comparable to U, then the singlet subspace has a ground state 
predominantly characterized by 15*2) configuration, which is well separated from the four 
excited states characterized by l^i), l^s), 15*4), \ Se). The doubly occupied IS'4) and the singly 
occupied l^i) configurations, which are connected by tunneling between dot 1 and dot 2, 
get very close in energy. When these states are degenerate, two eigenstates can be obtained 
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by \Ut,) = ^(l^i) ± 1^4)). Similarly, we get \Ut,) = ^(l^g) ± l^e)) from l^g) and \S,). By 
second order perturbation theory, the well-isolated ground state with dominant contribution 
from IS'2), has energy 

Ems = V - ( ^ + -] . (9) 

\A + V' + U -V + 2Vi A + V' + U-V-2ViJ 

We note that the singlet-triplet splitting is Ej^jt — Ej^js > for all range of A„ and t, and 
we have singlet as the ground state. 

Next, we consider the three electron states. In the fully spin polarized subspace, Sy = 
+3/2, there is only one state \S = 3/2, Sy = +3/2) = Cg^Cj^cj^lO), with energy given by 
i<;3/2 _ 2^ _|_ ^ _|_ 2V' + V. This state is characterized by having a spin up electron in each 
dot. For Sy = +1/2 subspace, it is composed of 9 singly and doubly occupied configurations. 
To simplify the qualitative analysis, we focus on a truncated basis composed of the following 
3 singly occupied configurations: \a) = Cg^Cg-j-c^lO), \b) = cl^cl^c\^\0) , |c) = Cg^Cg^cj^ |0), and 
two doubly occupied configurations \d) = Cg^cJ^c^ |0), |e) = Cgi-Cg^cj^ |0). Figure [2] shows 
resonance between configuration \b) and |e) when |A| = 0{U). The three-dimensional 
subspace with singly occupied configurations with 5*^=1/2 can be further decomposed by 
the total spin S, since 5* is also a good quantum number. For the subspace with S'=l/2, we 
use the Jacobi basis states Lq and Li:— 122 \Lq) = — |c)), \Li) = — 2\b) + |c)). 

For \Lq), the spin state in dot 1 and dot 3 is a singlet. For \Li), the spin state in dots 1 
and dot 3 can be written as a linear combination of triplets with Sy = and 5*^^ = 1. The 
remaining Jacobi state IL2) = + 1^) + |c)) is a total spin 3/2 state and is decoupled 

from all other states. In a similar fashion, we form Jacobi coordinates for the two doubly 
occupied configurations \X) = + |e)), and \Y) = -^(M) — In the subspace of 

S'=l/2 and 5'j;=l/2, with basis {l-^^o), \Li), \Y)}, the 3-electron Hamiltonian, 
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A + 2V +V -t 

-t U + 2V 

A + 2V' + V Vst 

V^t U + 2V 



(10) 



separates into the pair of Hamiltonians describing Jacobi basis states |Lo) and \Li) en- 
tangled with the doubly occupied configurations. Each sub-matrix can be diagonahzed 
and the eigenstates read: \Lq) = cos(0)|Lo) + sin(0)|X), \Lq) = sin(0)|Lo) — cos(0)|X), 
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|L+) = cos(^)|Li) + sin(^)|y), and = sm(^)|Li) - cos(^)|y), where tan(2(/)) = t/^, 
tan(2e) = VSt/^, and ^ = (A + 2^' - f/ - \/)/2. Here, we observe that each of the two 
Jacobi states, characterizing the (1, 1, 1) configuration, hybridizes with both doubly occu- 
pied configurations |X) and \Y) to form the eigenstates of a central-dot biased system. All 
four eigenstates, |Lq ) and \Lf), are current-conducting because electrons can be removed 
from the orbitals in the edge dots to make a transition from the three-electron state to a 
two electron (1,0,1) configuration. 

In Fig. [3](a) we show the evolution of the five lowest energy levels of the three electron 
complex in the Sy = 1/2 subspace as a function of bias A in the central dot. At A = 0, the 
spectrum is divided into 2 bands. The lower band consists of |Li), \Lq), and IL2) states, 
which are all characterized by singly occupied configurations. The upper band consists 
of states with dominant configurations \d) and |e). As A increases, the energy difference 
between the singly occupied configurations and specific doubly occupied configurations \d) 
and |e) diminishes. However, the ground state is always the \Lf) state in the figure. The 
blue curve represents the spin-3/2 state which does not interact with all other levels due 
to the conservation of total spin of the Hamiltonian. In the plot, the levels are artificially 
shifted for better visualization. In the inset of Fig. [3]^a), the proper energy levels around 
the anti-crossing point are shown in detail. Figure [3](b) shows the configuration content of 
the ground state as a function of bias A. At A = the ground state is dominated by singly 
occupied configuration but at higher bias, A ^ U, the doubly occupied configuration 
\Y) reaches around 50% content of the ground state. 



B. Current through a linear triple quantum dot 



Theory of sequential tunneling through a triangular TQD molecule has been described in 



detail in Ref. 



171 . Here we extend the approach to include both electron-phonon interaction 



and cotunneling and apply this theory to describe current and spin blockade in a LTQD. 
Following Ref. [iTI, current between lead r and a TQD device in the vicinity of a QP involving 
N = 2 and + 1 = 3 electrons can be written as a difference between the current from the 
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lead to the TQD and a current from the TQD back to the lead r: 

ir,cr aM,l3N+i 

+ E W^r'(/3iv+i^«Jv)^/3^+i, (11) 

where la^) is an A^-electron many-body eigenstate of the isolated TQD with energy E^j^ 
and associated steady state probability PdMi which is obtained by solving the rate equation, 
which is explained below. The sequential tunneling rate, W^'^'^{aN Pn+i), provides the 
rate of transition for the TQD from an N-electron state to an (A^ + l)-electron state 
due to first order perturbation from the lead r. Details of sequential tunneling rates will be 
provided later. 

The probabilities Pq-j^ '^ diagonal matrix elements of the reduced density matrix 

p. The time evolution of these diagonal matrix elements is described by the Pauli master 
equation, 

3 

Po.^ = E Y^P^N'^^f^^' ^ - Po^N^iaN (12) 

N'=2 Pi^, 

where transition rates Waf^^p^, are calculated using Fermi's Golden Rule. We consider 
sequential tunneling rate W^'^'^ in first order in coupling to the lead r, intra TQD phonon- 
induced relaxation rate W^^, and second order cotunneling rate W!^°^. The master equation 
is solved to obtain steady-state solution for the probabilities, Pa^i by setting the time 
derivatives to be zero. 

With the coupling to a lead r in Eq. ([2]), the first order sequential tunneling rates read 

2 

5 {uJali - erk) friy^af}), (13a) 

2 

(5 - erfc) (1 - /r(c^a/3)), (13b) 

where /r(e) = 1/ (exp[(e — fir)/kBT] + 1) is the Fermi function of the lead r, Uai3 = -E'/3^r+i ~ 
i?Q,^, and erk is the energy of a state associated with wave vector k of lead r. We remark that 
the summation over index i in the sequential tunneling rate refers to summing the tunneling 
contributions from the S and P orbitals in a quantum dot. By expanding the norms of the 
complex-valued matrix elements in above equations and introducing an integration variable 



seq 



2tx 



seq 



2tx 



(/^TV+l «Af) = Y 5Z 



^t[(A;)Q^|/3Ar+i) 
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w, the sequential tunneling rates can be also expressed as follows, 



N+l) 



2n ^ 



(14a) 
(14b) 



with 



generalized spectral functions^^ of 



the 



TQD, 



Xlo-('^Jv|cj(T|/3iv+i)(/3Af+i|cjo-|ttAf)'^('^ — Wq,^), and generalized spectral function of the 
lead r, B\-{uj) = J2k^i(^)(^]j(^))*^(^ ~ ^rk)- By substituting the sequential tunneling rates 
in Eq. flTTl) with Eq. flT^ . one can relate the current through a TQD with the spectral 
functions of the TQD and the leads. 

We now provide relaxation rates due to electron-phonon interaction and cotunneling. For 
large source-drain bias voltage \eVsd\ ^ \tij\, the change in the on-site energy of dots due 
to the source-drain bias will take the system off the resonance, away from the QP. In this 
regime, the current is dominated by inelastic tunneling between orbitals of neighbouring 
quantum dots due to electron-phonon interaction. The phonon emission-induced relaxation 
rate^ reads, 

2 



N 



2tx 



5 {Eaj^ - Ep^ - hwq) g{nujq, T), 

(15) 

where hwq = hcs\q\ is phonon energy, and gihujc^) is the thermal occupation number for 
phonon mode q at temperature T. Spin blockade occurs when the spin-3/2 polarized states 
I tta) become a trap state, with W,^^'^{a'i — )■ (32) = 0. However, the spin blockade can be 
lifted if we allow cotunneling. We consider cotunneling transition rate,— which involves 
an exchange of electrons between a lead r and the TQD in a spin 3/2 state, 

2n 



w: 



cot 



<j,(j ,k,k' 



a 



Q372 



Fas Eryr, £{,/ 



72 



(16) 



where Qg,^^ 



(Q;3|c|g.|72), I72) is a triplet state, and is the energy for an electron with 
wave vector k and spin a of the lead r. 
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III. TRANSPORT AND SPIN BLOCKADE 



In this section, we compute and discuss the transport properties and spin blockade in a 
LTQD at both SQP and AQP. We set on-site Coulomb repulsion between S orbitals to be 
Uii — U22 — Uss — U — 3.0 meV, and we use U as the unit of energy scale. For S and P 
orbitals in the same dot, we set U14 = U35 = U' = 0.94C/, and U44 = = U" = 0.96?7 
for P orbitals in the same dot. Wc set tij = t = —6.0 ■ 10~^ U for tunneling between S 
orbitals in neighbouring dots. We set tij — t' — —6.2 • 10~^ U for tunneling between the S 
and P orbitals on neighbouring dots. We set — V — 0.2 U between neighbouring dots 
and Vij — V — 0.1 U between dots 1 and 3. The energy difference between S and P energy 
levels. Asp, in the same dot is taken to be 0.8 U and 0.25 U in different cases considered 
below. 

The tunnel coupling for the tight binding chain in the leads is taken as ti = = —2.0U. 
The large tunnel coupling for the leads allows a wide energy band, which increases the 
amount of available states for transport. As for the dot- lead tunnel coupling t^, we set 

= -1.0 • and = -1.1 • lO'^U. Only the S and P orbital in dot 1 is connected 

to the left lead. Symmetrically, we set — if, and tf — t^. The rest of the tunnel 
coupling parameters are zero in our model. For interaction between electrons in the TQD 
and bulk LA phonons, we use the following GaAs parameters: A(g) = ^y D^hWq/2pc^, where 
D = 2.9U, p = 5300 kg /m^, Cg = 3700m/s and Ug = CgQ. 

We measure current in unit oi Iq — eltil"^ /h\tL\- We assume total potential difference eVgd 
across the two leads and a hnear decrease of this potential across the device. The chemical 
potentials on the two leads are given by /il — eVsd/'^ ^-nd pu — —eVsd/2. The on-site 
energies are given by £^1,(4) (V^^^) = + Ei.^ = E° , and E^^i^^^iVsd) = ^3 (5) -eV^^/e 

respectively, and electron temperature in all calculations is set to /c^T = 2.0 • 10~^ U. 

A. Quantum Interference-Based Spin Blockade 

We first consider transport through the SQP: {(1, 0, 1), (2, 0, 1), (1, 1, 1), (1, 0, 2)}, and 
we put P orbitals {i — 4,5) in dot 1 and 3, respectively. We set — E^ — —U — V 
and E2 — —2V' in order to bring the four charge configurations into resonance. For the 
present case, we set a high single particle level spacing Asp — O.SU in the edge dots. A 
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large energy spacing between the S and P orbitals allows one to focus on a few lowest states 
for transport at bias \eVsd\ ^ U. For instance, Fig. H] shows the energy diagrams of the 
relevant 2-electron and 3-electron states near the SQP as a function of Vgd- In the presence 
of small Vsd, the energy spectrum does not alter much and the wavefunctions remain similar 
to the wavefunctions at zero Vgd- The inset in Fig. H] categorizes the states associated with 
the energies in the main figure. There are four active 2-electron states: one singlet | S), and 
triply degenerate triplet states, |T^'°). These four states are characterized predominantly 
by the (1,0,1) charge configurations as discussed in Sec JII A[ For = 3 subspace, there 
are four spin-1/2 states below the spin-3/2 states. The four spin-1/2 states are \Lf) and 
\Lq) and their counterpart in the Sy = —1/2 subspace. In the absence of magnetic field, 
these states remain degenerate. Next up in the three-electron subspace are the quadruply 
degenerate spin-3/2 states, | S = 3/2, Sy = ±1/2, ±3/2). The last four levels are \L^) 
and |Lq) states and their counterpart in Sy = —1/2 subspace. We emphasize that these 
three electron states are admixtures of (2,0,1), (1,1,1), and (1,0,2) configurations with 
comparable weights except the spin-polarized states as discussed in Sec JII A[ Based on the 
analysis of wavefunctions, obtained from the exact diagonalization of a single-band Hubbard 
Hamiltonian, the only dark channels in the LTQD are the spin-3/2 states. As the spin-3/2 
wavefunctions, | S = 3/2, Sy = ±1/2, ±3/2), do not overlap significantly with the two- 
electron triplet states, | T^'°), when an electron is added or removed from the edge dots, the 
conventional spin blockade is not expected in this regime. 

Figure [5t^a) shows the current /(Kd) of a LTQD and Fig. [5](b) shows the steady state 
occupation probability of the four spin-3/2 states as functions of Vsd- This was done without 
the cotunneling effect. The / — Vgd curve is symmetrical with respect to the bias direction 
as it should be at SQP. The most prominent feature is that the vanishing of the current and 
therefore significant negative differential conductance associated with high occupation prob- 
ability of the spin-3/2 states. As shown in Fig. [5](a), the current is completely suppressed 
at a very limited bias regime, this is very different from the I — V curve in the spin blockade 
regime in a DQD. These numerical results are obtained from a five- level Hubbard model, 
and the negative differential conductance is not reproduced when we use just the three- level 
Hubbard model for the transport calculation. This implies that this negative differential 
conductance is related to the existence of the high-energy P orbitals. 

In order to explain this negative differential conductance, we need to study Sy = 3/2 sub- 
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space with all five orbitals. There are 10 possible configurations for three spin-up electrons 
in five orbitals. Using the Hubbard model with these five orbitals, the configuration with 
the lowest energy is |a) = c|^C2|C3||0), and the next two configurations are \ h) = c\^cl^cl^\0) 
and I c) = cj^ 63^04^10). Configuration | a) is separated from | b) and | c) by an energy gap of 
~ U +Asp — A. The other 7 configurations are even further away in energy. The Hamiltonian 
of this low energy configuration subspace in the basis of {| 6), | a), | c)} is, 

El + 2E3 + Asp + U' + 2V -t' 

-t' E1 + E2 + E3 + 2V' + V -t' 

-t' 2Ei + + + U' + 2V 

(17) 

where Ei and E-^ are almost identical when Vsd is small. This Hamiltonian matrix looks 
similar to the 2-electron triplet Hamiltonian, Eq.(l5]), except that the tunneling matrix el- 
ements acquire a negative sign for the three electron system. This negative sign is simply 
due to the anticommutation relation between fermionic operators. Exact diagonalization of 
the above Hamiltonian gives a ground state, 

\m = v(\a) + -\b) + -\c)], (18) 



where coefficients 771(2) are of the same order of magnitude as the coefficients 71(2) for |T+) 
in Eq.(l7]). This can be understood by analyzing the Hamiltonians. The energy difference 
between the configuration | Ti) and | T3) is given by | A + V — V|, and the energy difference 
between the configurations | b) and | d) is given by |A — Asp — U' — V + 2V'\. Considering 
that A and Asp are both of the order of U, the two energy gaps are actually comparable. 
In general, hybridization of configurations | i) and | j) in a wavefunction can be estimated 
by ^^te'- III our case, the S-P tunnel coupling t' is of the same order of magnitude as 
the S-S tunnel coupling t. This explains why 771(2) are comparable to 71(2) in magnitude. 
Furthermore, 771(2) have opposite signs with respect to 71(2) because the off-diagonal matrix 
elements in Eq.f ll7p and Eq.([S]) have opposite signs. Figure EI^c) presents the norm of 71 and 
rji from the exact diagonalization of the five-level Hubbard model as a function of Vsd- 

Next, we look at the rate equation for the state |3/2) when the system is subject to a 
positive source-drain bias, i.e., charging electron from left dot and removing electron from 
right dot: 

= -Wp{3/2 T+)P3/2 + WP{T+ ^ 3/2)Pt+ . (19) 
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Note that the only allowed 2-electron state is IT"*") because the total spin cannot change 
by more than 1/2 by adding an electron. The phonon relaxation does not play a role here 
because |3/2) and |T+) are the lowest energy states in their own spin- resolved subspaces, 
respectively. For simplicity, we have ignored the cotunneling contribution in this analysis, 
and numerical results in Fig. Oare also obtained without the cotunneling terms. Cotunneling 
effects will be discussed below. In order for 1 3/2) to be a trap state, the outgoing part of the 
rate equation should be almost equal to zero. The outgoing sequential rate is, approximately, 

WP{?,I2 ^ T+) = |(T+|tf(A:)c3t|3/2) + (T+|tf (A;)c5t|3/2) |^ 

k 

X 5(^3/2 - Et+ - ekR){l - fn), 

2 

6{Ey2-ET+-ekR){l-fR), (20) 

where if'{k) = e~*'^°/v^27r and fji is the fermi function for the right lead. The coefficients 
7i and rji are defined in Eq.Q and Eq. (fT8l) . respectively. The expression (|7i|tf — 1^11^5^) 
gives the interference between the two possible paths of removing an electron (via the S and 
P orbital) from the right dot. The minus sign in the expression stems from the fact that 
rji and 71 have opposite signs, and the origin of this sign difference was already explained 
immediately following Eq.(IT71). We see that the condition for the quenching of the sequential 
tunneling rate is |7i/'/7i| = /^f |- Figure[5]^d) presents the ratio |7i/77i| as a function of Vsd- 
At points of strongest current suppression, we observe that the ratio indeed matches the ratio 
of \t§/t^\- In short, the negative differential conductance sets in whenever the two possible 
paths of electronic transport become comparable in amplitude and interfere destructively. 
This destructive interference is possible only for the transport channels through spin-3/2 
states. In terms of spin configurations, the transport channel |3/2) — )■ |T+) involves the 
two paths (ti,t2,t3) (ti5t2) and (ti,t35t5) i'tiA-i), which can destructively interfere. 
For all other transport channels, electronic transport occurs with much higher probability 
amplitude via the 5* orbital in the edge dots at low source-drain bias. The existence of the 
dark channel through |3/2) makes the TQD molecule to be trapped in |3/2) state. 

Figures [n](a) and (b) show the current through the LTQD and the steady state probability 
distribution for the spin-3/2 states in the parameter space of {Ei = E^, E2) at a small bias, 
eVsd = 0.01 U, respectively . In this calculation, the cotunneling effect is included. Although 
the quantum interference-based spin blockade is formed under a very specific condition. Fig. 

16 



2n y--^ 



\ii\tz-\ni\t^) 



-ika 



Mjo) shows that the interference-based spin blockade can still be observed in the parameter 
space of on-site energies. The cotunneling effects can be analyzed when we add terms 
E.=L,i?E& Wr\(33 ^ 3/2)P/^3 _ Wr\S/2 ^ /33)P=^/^ where | /^s) is a three- 

electron state with S = 1/2, in Eq. flTU]) . In Eq.f l20p . we analyze the condition for the 
transition rate from | 3/2) to | T+) to vanish. With cotunneling included in the model, we 
should analyze the condition for the transition rate from | 3/2) to each | P^) state to vanish. 
In principle, each transport channel has its unique condition for the quenching, and the 
interference-based quantum spin blockade will be lifted. However, the additional rates due 
to cotunneling are much smaller in amplitude as they scale with Itfl^ for the second order 
processes. The system still gets blockaded in the spin-3/2 state when the sequential tunneling 
driven transition {W^'^{3/2 — )■ T"*")) vanishes, because the incoming rate iy^'^'^(T+ — )■ 3/2), 
a first order process, scales with |tf p and is around 5 orders of magnitude larger than the 
rates driven by cotunneling processes. 

B. Symmetrical and Asymmetrical Spin Blockade 

Next, we look at the SQP again with a different single-particle level spacing. Asp = 0.25 
U. In this case, we will consider a wider range of source-drain bias with eVgd > U . Figure [Tl^a) 
shows the current, /(Kd); of the TQD near the SQP. We again have a symmetric liysd) with 
respect to the bias directions and, therefore, the observed negative differential conductance 
is also bi-directional. We will focus on the positive bias direction for the following discussion. 
We note that there are 2 regions where the current is strongly suppressed in the positive bias 
direction in Fig. [Tl^a). One point is at the low bias regime, eVsd ^ U , and the other point 
is at the high bias regime such that on-site triplet occupation is allowed in the transport 
window. From Fig. [Tl^b), we see that the system is trapped in (1,1,1) spin 3/2 states 
whenever the current is significantly suppressed in Fig. [Tl^a). The strong current suppression 
at the low bias is due to the quantum interference-based spin blockade we described in the 
previous section. As source-drain bias is further increased, the wavefunction inside the LTQD 
also changes. Gradually one path of electronic transport becomes preferred and quantum 
interference vanishes. At the high bias, the second current suppression is identified to be the 
more familiar spin blockade phenomenon in the double quantum dot, and it is characterized 
by an extended region of current suppression over a wider range of source-drain bias. At 
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high bias, hybridization of levels becomes insignificant, and it is instructive to look at each 
eigenstate as a particular localized configuration. Figure [8] presents a schematics of how this 
high bias spin blockade is formed and lifted in the LTQD at high bias. The spin blockade 
is formed when the on-site triplet becomes accessible in the left dot but not in the right 
dot in the transport window when a positive bias is applied. Due to the phonon-induced 
relaxation, the on-site triplet in the left dot will relax by allowing electron-phonon scattering 
to re-distribute the electron from the P orbital in the edge dot onto the 5* orbital in the 
central dot. When the on-site triplet state in the right dot is still too high in energy for 
occupation, the system gets stuck in this (1, 1, 1) spin-3/2 configuration. This spin blockade 
is lifted when the bias is further increased so the on-site triplet become accessible in the right 
dot too. Then the phonon-induced relaxation will again help transfer the electron from the 
central dot onto the right dot. We remark that the spin blockade does not happen in this 
model if the phonon-induced relaxation mechanism is removed. From this picture, we can 
derive the spin blockade regime from the parameters we used. The energies of the relevant 
configurations are 

^(tits) = E^ + E^ + V = -2.1U , 
E(tit3t4) = 2E, + + ^3 + f/' + 2y + ^ = -imU + ^ , 

D D 

^(tit2t3) = E^ + E2 + E^ + 2V + V' = -2.1U , 

^(titsts) = Ei + 2Es + Asp + U' + 2V-^ = -IMU - ^ . 

o o 

For an electron to move from the left lead to the TQD, £'(^1^3) + ee=-E(tit3t4) fo^' 
electron energy Se < Hl = Thus, we get eVsd > 0.57t/. -E(tit2t3) is always lower 

than E(tit3t4) for forward bias, so the relaxation from | tit3t4) to | tit2t3) is allowed. For 
the spin blockade to occur, the transition from | tit2t3) to | tit3t4) should not be possible 
by phonon emission. So, i?(tit2t3) < -^(tit3t5)5 which leads to eVsd < 1.14f/. Therefore, 
the spin blockade regime is 0.57t/ < eVgd < 1.14t/, which agrees very well with the numerical 
result in Fig. [7](a). 

Next we consider the current of the LTQD at the AQP: (Oil), (012), (021), and (111). We 
again use a five-level Hubbard Hamiltonian for transport calculation. We put S orbital in the 
left dot, and S and P orbitals in the central dot and the right dot. In weak tunnel coupling 
limit, the 4 charge configurations should be on resonance, and we set Ei = —V — V and 
E2 = E3 = -U - V. Figure EKa) shows the /(Kd) of the TQD near the AQP. As expected, 
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the current of the TQD near the AQP is very different under the two bias directions. Figure 
Mjo) confirms the spin blockade where the current is severely suppressed in the positive bias 
direction in Fig. IHl^a). This phenomenon is in close analogy to the case of a DQD, and can 
be easily explained. In the positive bias, electron is injected from the left dot. The transition 
from (0, 1, 1) triplet to a (1, 1, 1) spin 3/2 state does not require the formation of on-site 
triplet. In the negative bias direction, electron is injected from the right and transition from 
(0, 1, 1) triplet state to a spin 3/2 state requires the formation of an on-site triplet in the right 
dot. Therefore, before the bias threshold = E{0, 1,2*) — E{0, 1, 1), no spin blockade is 
expected to be formed. (0, 1, 2*) represents charge configuration in which one of the electrons 
occupies the P orbital in dot 3. When the applied bias exceeds the threshold, the condition 
E{0, 1, 2*) > E{0, 2*, 1) > E{1, 1, 1) is also satisfied. Thus, either due to resonant tunneling 
or inelastic process, this additional electron can always be removed from the right lead. 
Therefore, there is no spin blockade in the negative bias direction. We remark that the spin 
blockade at AQP can be formed without the assistance of any relaxation mechanism. So 
the spin blockade of a TQD at AQP is almost identical to the spin blockade in a DQD. 
Figure fTOl presents the transport triangle in the parameter space {Ei, E2 = E^) at a positive 
bias. In this figure, the light trail at the tip of the triangle is proportional to the on-site 
singlet-triplet gap in the central dot. This transport triangle, although generated under 
the specific condition E2 = E-^, provides similar information that one can extract from the 
transport triangle for the DQD. 



IV. CONCLUSION 



We presented a theory of electronic properties and transport through a LTQD around 
QPs. We showed that the spin blockade could serve as a spectroscopic tool for the detec- 
tion of different spin states. Two different QPs containing the (1, 1, 1) configuration were 
discussed. A multi-band Hubbard model with five levels was used to describe the electronic 
properties and investigate the spin blockade phenomenon in the LTQD. For both QPs, strong 
current suppression and negative differential conductance were predicted. At the SQP, sup- 
pression in conductance was obtained under two different source-drain bias regimes. When 
the bias is small, the electronic transport involving spin-3/2 states takes place either via 
the 5 or P orbitals in the edge dot with comparable amplitude and results in a destructive 
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interference. In high bias regime where electron tunnels onto the P orbital in the edge dot, 
spin blockade is facilitated by spin-conserving relaxation mechanisms, such as interaction 
with LA phonons studied here, and formation of the trap state. At the SQP, the spin block- 
ade phenomenon is bi-directional, in contrast with the spin blockade in a DQD. We also 
discussed spin blockade at the AQP. The spin blockade formation and lifting in this case 
is in close analogy to the DQD case. The formation of the spin blockade does not involve 
any on-site triplets, only the lifting of the spin blockade requires the access to the on-site 
triplet states in the transport window. Similar to the DQD, the spin blockade phenomenon 
at AQP only occurs only in one of the bias directions. 
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FIG. 1: (Color Online) (a) Schematic picture of a LTQD with one electron spin each connected to 
leads. The leads are modelled with ID tight binding chains, (b) The TQD in (1, 0, 1), two-electron 
configuration, (c) Schematic picture of single particle energy spectrum of a TQD when the central 
dot is biased with A. The gap A^p denotes the energy difference between S and P orbitals on a 
dot. 
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(a) 




FIG. 2: (Color Online) Two resonant three-electron configurations in a centrally biased TQD : 
(a) one of the (1,1,1) singly occupied configurations and (b) one of the (1,0,2) doubly occupied 
configurations. 
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FIG. 3: (Color Online) (a) The three electron molecule energy spectrum with total Sy = 1/2 as 
a function of bias A, obtained from the single-band Hubbard model. As A increases, the five 
energy levels anti-cross. The blue curve, corresponding to total spin-3/2 state, does not interact 
with the other states. The energy levels are artificially shifted by a constant values for better 
visibility. The inset shows the energy levels near the anti-crossing point, (b) shows the projection 
of the ground state onto \Li) (the red curve) and \Y) (the black curve) states. The wavefunction 
of the ground state is obtained from exact diagonalization of the single-band Hubbard model. The 
blue-dashed lines provide the same information but obtained from the analytical approximation 
for \L'^) = cos{6) I Li) + sin(0) | Y) discussed in the text. 
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FIG. 4: (Color Online) The energy spectrum of a multi-band Hubbard model as a function of 
source-drain bias. In the figure, the blue curve represents the four 2-electron states. They are very 
close in energy and looks degenerate in the energy resolution present here. The green curve is the 
quadruply degenerate spin-3/2 states. The red curves are each doubly degenerate spin-1/2 states. 
Inset: A summary of the states involved in the main figure at a particular value of Vgd- 
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FIG. 5: (Color Online) (a) The current KVsd) as function of the applied source-drain bias Vgd of a 
LTQD at the SQP in the low source-drain bias regime. Note zero current at Vgd ~ ±0.08. (b) The 
steady-state occupation probability of the spin-3/2 states as a function of Vgd- Panels (a) and (b) 
together indicate that the spin-3/2 states are related to the bi-directional, quantum interference- 
based dark channel in a LTQD. (c) Projection of the triplet state |T+) onto the configuration 
c|[|C2^|0), and the projection of the spin-3/2 state |3/2) onto the configuration cJ^Cg^Cg^jO). (d) 
Ratio of matrix elements | ^ | (see text for the definition) in the unit of the ratio | ^ | . Spin blockade 
is formed when the red curve intercepts y = 1 line in the figure. Panels (c) and (d) are presented 
to illustrate the formation of the spin blockade in the positive bias direction. 
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FIG. 6: (Color Online) (a) Current of the LTQD in the parameter space of {E\ = £'3,i?2)- The 
transport region manifests a rounded boundary, which indicate the states involves in the electronic 
transport are highly hybridized states. The transport region is separated into two parts by a 
thin line of strong current suppression. This is the region of dark channels, (b) The steady-state 
occupation probability for spin-3/2 states. High amount of spin-3/2 states are found exactly where 
the currents vanishes in Panel (a). 
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FIG. 7: (Color Online) (a) Current of the LTQD at SQP as a function of Vgd- At high bias (in 
both bias direction), we observe a robust negative differential conductance, (b) The steady-state 
occupation probability for spin-3/2 states. The peaks near the low- the figure Indeed, the current 
suppression is associated with the spin-3/2 states. The much more extended current suppression 
is due to the spin blockade phenomenon. 
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FIG. 8: (Color Online) (Schematic representation of lifting of spin blockade, (a) (2,0,1) configura- 
tion obtained from the (1,0,1) configuration when an additional electron tunnels onto the P orbital 
of dot 1. (b) Due to phonon- induced relaxation in the model, the added spin moves from the P 
orbital in dot 1 to the S orbital in dot 2. However, it does not proceed further to the dot 3 because 
this costs energy, (c) At larger source-drain bias in the positive direction, the energy levels in dot 
3 are lowered with respect to that of dot 2. Thus, phonon-induced relaxation assist the electron to 
move onto dot 2. 
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FIG. 9: (Color Online) (a) The current of the LTQD at AQP as a function of Vgd- The current 
response is asymmetrical with respect to the bias direction. Similar to a DQD, current suppression 
is only observed in one direction of the bias, (b) The steady-state occupation probability for 
spin-3/2 states. The current suppression is associated with the spin-3/2 states. 
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FIG. 10: (Color Online) The current of the LTQD at AQP in the parameter space {Ei, E2 = E3). 
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